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Abstract 

We study the (2+l)-dimensional single-step model (SSM) with a tunable parameter p. Using 
extensive numerical simulations, we find dependence of interface exponents and scaling exponents 
derived from iso-height clusters on p. Also we test the iso-height contours for being Schramm- 
Loewner evolution (SLE K ) curves. We find that there is convincing evidence that they are with 
k ~ 8/3 at p = 0, but for p ^ 0, the results do not match with each other. Furthermore, the 
boundary conditions become important for p / 0. The right choice of boundary conditions at 
p = 0.5 such that Gaussian free field universality class emerges, remaining unknown. 

I. INTRODUCTION where = means statistically the same, and 

a > 0. A possible way for classification of 
Growth processes and rough surfaces are wious gurface grQwth modelg ig baged Qn 

interesting topics in physics from both theo- scaling behaviour of surface width> w ^ L) = 



ELQ. 



retical and experimental points of view Ufciaj- ^([h^t) - (h)} 2 ), and (••) denotes spacial 

Various discrete models have been suggested averaging For a non . equ ilibrium growth, the 

to describe surface growth processes, for ex- width ig ted to have the following scal . 

amples see |2|. These models produce a j ng f orm 
self-affine interface h(x) such that its proba- 
bility distribution function remains invariant 

w 2 (t,L)~L 2a f(-). (2) 

The scaling function / usually has the asymp- 
totic form f\x — > oo) = constant and f(x — > 
0) ~ x 213 . The saturation time t s has the scal- 
ing ansatz t s ~ L z . The universality class can 



under scale transformation: 



/i(x) = b- a h(b*) 
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and growth exponents, a and p, respectively. 



The dynamic exponent is z = ot/fl. A new tool for study of domain walls in 

For some equilibrium-rough interfaces, the critical systems is the theory of Schramm- 
width behaves logarithmically Q, i.e., Loewner Evolution (SLE) {9]; for a review 



w 2 (t,L) ~ hit for t <C t sat , and w 2 (t,L) ~ see [lOj. The scaling behaviour of the twa- 
in L for t 3> t s , with z = 2. dimensional critical lattice models can be re- 
Whether rough surfaces fall into dis- fleeted in the statistics of non-crossing ran- 
crete universality classes is unknown, how- dom curves which form the boundaries of 
ever many self-affine surfaces are observed clusters on the lattice. 

to belong to certain universality classes. In the 1920s, Loewner studied simple 
Two well known universality classes are de- curves growing from the origin into the up- 
scribed by continuous Langevin equations per half-plane EI [11] . Loewner's idea was 
i.e., the Edwards- Wilkinson (EW) {4] and to describe the evolution of these curves in 
the Kardar-Parisi- Zhang (KPZ) |j] equa- terms of the evolution of the analytic func- 
tions. The KPZ equation is given by tion gt{z), which conformally maps the region 
gur-x +\ \ H\iT t to H, where K t is the hull (the union of 
dt = " V ' h + 2 I™!' + '<»•«>• < 3 > the curve and the set of points whieh can not 
where the relaxation term is caused by a sur- be reached from infinity without intersecting 
face tension v, and the nonlinear term is due the curve )- He showed that this function sat- 
to the lateral growth. The noise rj is uncorre- isfies the following differential equation: 
lated Gaussian white noise in both space and dg t 2 
time with zero average i.e., (?y(x, t)) = and 9t\Z) — it 
(i]('x,t)rj{x' ,£')) = 2D5 d {yi — x.')5(t — if). For for a real continuous function £ f , related to 
A = 0, the EW equation is retrieved whose the image of the tip of the curve under g t . 
exact solution a = (2 — d)/2 and z = 2, is More recently, Schramm proposed the idea 



known in (d + l)-dimensions [4|. For KPZ that a measure on the continuous driving 

equation, due to additional scaling relation functions & would induce a measure on the 

a + z = 2, there remains only one indepen- set of growing curves in H, and showed that 

dent exponent, say a. The exact solution the latter measure is conformally invariant 

only exists in Id [5[ which gives a = 1/3. if and only if £ t = y/RB t . The diffusivity 

In 2d, the exponent is available only by var- constant k determines the critical exponents, 

ious simulations and theoretical approxima- hence, the universality class of the system 

tions ranging from a = 0.37 to 0.4 p-|8(. in question. Some authors have argued that 



SLE K may be app 



nomena as well 



12 



led to surface growth phe- 
-|l5|. 



3er we study 2d single step 



In this pa 
model (SSM) 
(SOS) model 



16Hl9| , a kind of solid on solid 
20J which is defined as follows. 
One column is chosen randomly, if it is 
a local minimum then its height is increased 
by 2 with probability p + ; one can consider 
this process as deposition. If it is a local 
maximum then its height is decreased by 2 
with probability p_, as a process of desorp- 
tion or evaporation. This definition guar- 
antees that in each step, the height differ- 
ence between two neighboring sites is ex- 
actly one. Overhanging is not allowed in 
this model and the interface will not develop 
large slopes. We start with the initial con- 
dition h{i,j;t = 0) = [1 + (-l) i+J ]/2 where 
1 < i < L x and 1 < j < L y . We impose the 
condition p + + p_ = 1, so due to up/down 
symmetry p + <H- (1 — p + ), we just need to 
consider p + < 0.5. Here after we shall refer 
to p + as p. 

Several papers have investigated this 
model, and claimed that for p = 0.5 (p ^ 0.5) 
the model belongs to EW (KPZ) universal- 



ity class 



17|-ll9|- Plischke et al. 17], showed 



that for p = 0.5, this model in Id is reversible 
and can be exactly solved by mapping to the 
kinetic Ising model. They found a = 1/2 
and z — 2. Furthermore for p ^ 0.5, they 
mapped the interface model onto the driven 



hard-core lattice gas, and focused on the av- 
erage slope of the interface, so in an approx- 
imate way they showed that the equation of 
average slope is in agreement with the Burg- 
ers's equation [121]. Thus claiming that the 
universality class is that of KPZ equation for 
p 7^ 0.5. They also simulated this model 
for p = 0.0, 0.25, and claimed that in the 
limit large L, the exponent z converges to 



zkpz = 3/2. Simu 



ations by the same au- 



thors on SSM in 2d 18[] , provided the scaling 
exponents a ~ 0, z ~ 2 for p = 0.5, and 
a f=s 0.375, z ~ 1.64 for p = 0. Kondev 
et al. [19I ]. simulated SSM on the square lat- 
tice with size L = 128. They confirmed that 
p = 0.5 shows Gaussian behaviour whereas 
p — 0.1 shows KPZ behaviour. However they 
found that p = 0.3 consistently resembles 



17| . attributing 



p = 0.5 contrary to claims of 
this to the above mentioned crossover from 
initial Gaussian behavior. 

We ran large simulations in 2c? for SSM, 
and estimated interface exponents and scal- 
ing relation for the iso-height clusters and 
contours. Our findings are consistent with 
the above within numerical errors. Further- 
more we checked SLE evidence for the iso- 
height contours and found that for p = 0, the 
curves are described by SLE 8 / 3 . But for other 
values of p, the results of SLE evidence, are 
not compatible with each other. Our guess 
is that the iso-height curves can be described 



by SLE, if we choose correct boundary con- x on scale b, is 19] 
ditions for interface. For this purpose, we M 

have checked the SLE evidence for domain C b (x) = ^ K x + &e ™) ~ M x )l • ( 5 ) 

m=l 

walls of critical Ising model under two dif- 
ferent boundary conditions, and the results The °ff se t directions {e m }m=i are a fixed 
were consistent se ^ °^ vec t° rs summing to zero. In our 

™. . • i r 11 • case (square lattice), |e m | pointing along the 

lhis paper is organized as follows, m sec- ' L 3 

rm , ., , ,. , {10} type directions. For a self- affine surface, 

tion [11] we describe our simulations and the 

interface exponents, in section HTT1 we describe 

the ensembles of the iso-height contours, and 

in section JV] we test for SLE curves for sev- 



the curvature satisfies the relation 



19j: 



(C b (x) q ) ~ b a « with a g = qa, (6) 



eral values of p. where (■■) denotes spatial averaging. To 

check this relation, simulations were done on 
square lattice of size L = 10 3 with more 

II. INTERFACE EXPONENTS than 10 4 height configurations. We applied 

periodic boundary conditions in both direc- 

We ran simulations on a square lattice t ions. The results for three different values 

of size 50 < L < 700, averaging w(t, L) Q f g> are plotted in fig. ffl As expected from 

over more than 200 independent runs. Size self- affinity, the three curves should coincide 

L = 4000 was used only for computing 0. To with each other, which is the case within er- 

calculate the exponents of iso-height clusters ror bars. We have excluded the point p = 

and loops, more than 10 4 height configura- 0.5, q = 3, since curvature vanishes, 

tions were undertaken on square lattice with The interface exponents a, (3, z for various 

size L = 1000. Also in the part of SLE, strip values of p are p i ot ted in fig. E3 For p = 0, 

geometry L x = 3L y and L y = L was consid- these exponents are in good agreement with 

ered with 100 < L < 1000. We define one the KPZ universality class (g-S|. As expected 

time step as equivalent to L 2 tries for depo- t he values of a and [5 decrease with p, until 

sition or evaporation. p reacri es 0.5, where EW exponents are ob- 

To check the efficacy of our simulations, served. Also near p = 0.5, we have checked 

we calculated the roughness exponent a that the interface exponents satisfy the loga- 

from scale- dependent curvature in saturation rithmic scaling relation. So far the observed 

regime. The definition of curvature C&(x) at behaviour at p = and 0.5 is in agreement 



0.4- 



0.3- 



0.2- 



0.1- 



0.0 



1 1 

* = = -k - ; 


i i i i 




-ft 

\ -o- q = 2 . 




M _ . H _ q = 3 . 




\\ -1^-0 = 4 



















T Q 



0.1 



0.2 0.3 



0^ 



-r 

0.5 



FIG. 1: Three plots correspond to a q /q 
versus p. As expected from self- affinity, the 

three curves should coincide with each 
other, which is the case within error bars. 

with known results. Nevertheless, the behav- 
ior of these exponents in the intermediate re- 
gion is unclear (see fig. [2]). We observe that 
z increases to ~ 2.04 while p is increased to 
0.25. For p > 0.25, the exponent z remains 
constant at ~ 2. This is however not conclu- 
sive due to finite-size effects and small simu- 
lations. 



III. STATISTICS OF CLUSTERS AND 
LOOPS 

In this section we look at the iso-height 
contours and clusters which can give ex- 
tra information regarding self-affine inter- 
faces 
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21| • 



Consider an ensemble of height config- 
urations in the saturated regime and look 
at contours at a specific height say hs = 




FIG. 2: (Color online) Interface exponents 
a, (3 and z versus p. Simulations were ran on 
a square lattice of size 50 < L < 700. Size 
L = 4000 was used only for computing (3. 



(h) + Sy/([h(x) - (h)) 2 ), where 5 is a small 
real number setting how far around the aver- 
age height we are cutting the interface profile. 
Each island (or cluster height) is defined as 
a set of connected sites with positive height 
which identified by the Hoshen-Kopelman al- 
gorithm 22j. In this section we choose 5 = 0. 
The contour loops were uniquely determined 
by the algorithm explained in [23|]. A contin- 
uous geometric change is shown in fig. [3] by 
snapshots of the positive height clusters for 
various p. The ones with lower p have much 
dense clusters while for higher p the clusters 
become more porous and scattered. 

Fractal dimensions: One can expect (from 
self-similarity of clusters) a relation between 
the average mass of a cluster M and the ra- 
dius of gyration R, i.e., M ~ R Dc , where D c 
is the fractal dimension of clusters. The seal- 
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FIG. 3: (Color online) Snapshots of positive iso-height clusters for various p on a square 
lattice of size 200. The cut is made at the average height. 



ing relation between the average length of a 
cluster boundary I and it's radius of gyra- 
tion r, is given by I 



19 



21|. More- 



over, the relation between the average area 
of a loop and it's perimeter is I ~ a da where 
d a = df/2 (for clusters without holes). We 
took more than 10 4 height configurations on 
square lattice with size L = 1000 for averag- 
ing. Among them, we only present the plots 
of l(r) for several p in fig. HI The correspond- 
ing exponents are reported in fig. El We see 
that the scaling exponents such as df and D c 
for p = are similar to the ones for KPZ 



model 



13 



15]. When p is increased, the 



value of df (D c ) increases (decreases) until 
it approximately reaches 1.5 (1.91), in good 
agreement with corresponding exponents of 



EW model [13j, |24|. We have checked that 
the value of the exponents does not depend 
on the level of the height cut, although the 
range of scaling is affected after a certain dis- 
tance from a critical level height 25). In 
order to see the finite-size effects, we have 
also considered the reevaluation of the expo- 
nents by going to the lattices of larger size 
up to L = 3000 and by taking averages over 
a number of 2500 saturated height configura- 
tions and obtained similar results within er- 
ror bars. 

Distribution Exponents: We now look at 
the distribution functions of different statis- 
tical observables of the height clusters and 
contours, such as the contour length distri- 
bution n(l), cluster size distribution n(M) 




~ I — I I I I — 

10' 



-I 1 — I — I — I I I I 



10 2 



FIG. 4: (Color online) The average length of FIG - 6: ( Color online ) Distribution function 



a cluster loop boundary versus the average 
radius of gyration on square lattice with size 
L = 1000. Averages are taken over more 



of logr versus r (loop's gyration radius) on 
square lattice of size L = 1000, for various 
values of p. The slope gives the exponent 



than 10 4 height configurations, each of them T r ~ 1 shown in fi S- For more clarit y> the 



containing several iso-height clusters. 



curves are suitably shifted. 
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FIG. 5: (Color online) Scaling exponents of 
iso-height loops versus p. 

and distributions for the radius of gyration 
of the contours and clusters, i.e., n{r) and 
n(R), respectively. These distributions are 
expected to follow the scaling forms, n(l) ~ 
l~ T \ n(M) ~ M~™, n(r) ~ l~ Tr and n(R) ~ 



R~ TR 19|, |2l|. Simulations were run on 
square lattice of size L = 1000, to get more 
than 10 4 height configurations for 5 = 0. 
The distribution function of log r versus r is 
shown in fig. [6j 

All distribution exponents are reported in 
fig The dependence of these exponents on 
p is evident. The exponents r# and r r coin- 
cide within error bars. In the following, we 
will show that these distribution exponents 



depend on leve 



Olami et al. 



26) 



cut 5, previously noted by 



Dependence of exponents on the level of 
the cut: We now look at the dependence of 
our results on level cut 5. All previous re- 
sults were obtained at 5 = 0. We ran a num- 
ber of simulations for several values of 5, and 
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FIG. 7: (Color online) Distribution 
exponents (see text) versus several p. 
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FIG. 8: (Color online) The island-size 
distribution exponent as a function of the 
level cut for various p. 



2.8- 



concluded that the exponents such as frac- 
tal dimension of contours df and clusters D c , 
do not depend on the level cut. Also we com- 
puted the exponents for the island size distri- 
butions function of 6. We observe in fig. [8] 
that for p = 0.5 and p = 0.3 (and for p = 0.4 
as well not shown in the figure), the cluster- 
size exponent has a bowl-like functionality 
to 5. For large p, it has a much convexity 
around a certain level 5 m i n where the expo- 
nent reaches to its minimum value. However 
for p < 0.25, the exponents tm monotoni- 
cally decreases with 5. We have also com- 
puted Ti(S) (fig. [7] ) and T r (S) (not shown), 
and obtained the same behaviour. Observe 
that error bars are larger for 5 < 0, since in 
this region the number of positive iso-height 
clusters are scarce. 
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FIG. 9: (Color online) The exponent 77 
describing the power-law behaviour of the 
distribution of the cluster boundary length 
as a function of 5. 

IV. EVIDENCE FOR SLE CURVES 

In this section, the results are based on a 
series of extensive simulations of SSM on a 
strip geometry of size L x x L y with L x = 3L y 
and L y = L, in accord with the postula- 
tion of dipolar SLE. For each height con- 
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FIG. 10: (Color online) The length of 
spanning curves versus size of lattice 
L y = L; slope of each line (correspond to 
the each p) yields df(p). 

figuration, we find all spanning clusters at 
level 5 = in y direction, and assign cor- 
responding coastlines that connect the lower 
boundary to the upper one. We considered 
L = 100,200,300,400,500,1000 and exam- 
ined the scaling relation I ~ L f, in order to 
compute the fractal dimension of the span- 
ning curves. We could gather a number of 10 4 
sample curves from an approximately 7500 
independent saturated height profiles. The 
results of the fractal dimensions for various 
p, are summarized in fig. [10] and the first col- 
umn of table fl] They show a continuous de- 
pendence of df on p ranging from 4/3 to 3/2 
for p = to 0.5, respectively. For confor- 
mally invariant curves, the fractal dimension 



is related to t 



d f = l + k/8 <J Q. 



re diffusivity k by the relation 



27|, 



pute the winding angle 9 as defined by 
Wieland and Wilson |29] . For each curve 
we attribute an arbitrary winding angle to 
the first edge (that is set to be zero). The 
winding angle for the next edge is then de- 
fined as the sum of the winding angle of the 
present edge and the turning angle to the new 
edge measured in radians. The variance of 

the winding angle is believed to behave like 
I — I 

(6 2 ) ~ a + 61nL |29|, where for conformal 
curves b = | = 2(df — 1). We have computed 
the variance of winding angle of an ensem- 
ble of spanning iso-height curves for different 
p as a function of lattice-width size L, and 
confirmed that is linearly proportional to its 
logarithm with a universal coefficient b which 
depends on p (see fig. [11] and table [I]). The 
two fractal dimensions obtained from direct 
measurement and the one derived from the 
coefficient b, i.e., df = 6/2 + 1, are plotted in 
fig. [12] for a comparison. They almost coin- 
cide for p > 0.25 but with a small deviation 
for p < 0.25. 



To investigate if this relation holds, we corn- 



Let us now measure the diffusivity k from 
a direct SLE test by statistical analysis of the 
driving function £ t of the spanning curves. 
For an ensemble of proposed dipolar SLE 
curves with hull K t , in a strip geometry §a = 
{z 6 C,0 < < it A} of width 7rA, the evo- 
lution of the conformal map Qt{z) which maps 




FIG. 11: (Color online) The variance of 
winding angle versus logarithm of the lattice 
width for various p. 
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FIG. 12: (Color online) The fractal 
dimensions obtained from scaling relation 
/ ~ L d j- (squares) compared with the one 
derived from the slope of the linear plots in 
fig- ED (circles). 



^>A\K t to §a is given by [30 1 
dg t 1/A 



(7) 



dt tanh(((/ t (z) - &)/2A) ' 
with initial condition go{z) = z. In order 
to calculate sequences £(tj) for the iso-height 
curves, we used the algorithm introduced 



p 


dt K = 8(dt - 1) b k = 4b 


0.00 


1.339(4) 2.71(3) 0.70(1) 2.80(4) 


0.10 


1.332(4) 2.66(3) 0.69(1) 2.80(4) 


0.20 


1.350(9) 2.8(7) 0.74(1) 2.96(4) 


0.25 


1.415(7) 3.32(6) 0.82(1) 3.28(4) 


0.30 


1.465(6) 3.72(5) 0.92(1) 3.64(4) 


0.40 


1.495(5) 3.96(5) 0.996(4) 3.98(2) 


0.50 


1.501(4) 4.01(2) 1.01(1) 4.04(4) 



TABLE I: First column is the fractal 
dimension of spanning curves df and second 
column is k that is derived from df. The 

third column b is the coefficient of the 
variance of winding angle (6 2 ) plotted in 
fig- EH followed by the derived values of k 
from b in the last column. 



m [31J- First we consider the sequence of the 
points of a spanning curve as z® = (x°,?/j ), 
i — • • • N; the upper index represents time 
step, where t = and £ = £(£o) — 
0. At the ith step, we map the sequence 
{^ij, zl~ l , ■ ■ ■ z'n X } to the transformed and 



shortened sequence [z*, z\ +1 , • • • z l N }, by 
ing the map appropriate for dipolar SLE, 



us- 



2(L/7r) 2 log [cos(A,)], 



Z 3 



lf cosh[7r(^- 1 -g i )/2L] 



cos Ai 



h 



(8) 



j-i 



. For SLE 



where Aj = iryl /2L and ^ 
curves, the driving function satisfies £(t) 



tion. 



10 



3 OH , where B+ is the Brownian mo- 



10 



Another test of conformal invariance, is to 
compute the left passage probability P(p,<p), 
defined as the probability that a curve in 
the upper half-plane, passes to the left of a 
given point at polar coordinates (p, 0) . This 
is given by Schramm's formula [321 ]: 



+ 



^-COt(0) X 



^Q^fj-cot 2 ^)), 



(9) 



where 2F1 is the hypergeometric function. 
Although this equation is originally obtained 
for chordal SLE curves, it also holds for dipo- 
lar SLE with p < L [31]. 

We considered a substrate of strip geome- 
try with size L x /3 = L y = L = 1000, and ex- 
tracted all spanning contour lines at the level 
5 = 0, with total number of 10 4 . We then 
extracted the underlying driving function for 
each value of p = 0,0.1,0.2,0.25,0.3,0.4,0.5 
and examined if the result is in accord with 
the postulates of conformal invariance and 
that of predicted by Schramm's formula. We 
found that (see fig. [[3]) for p = 0, the diffu- 
sivity derived from driving function is k = 
2.56(6) and from left-passage k = 2.55(8) 
which is approximately compatible with frac- 
tal dimension of curves reported in table [H 
This is very close to the diffusivity of sel: 



avoiding walk (SAW) with k = 8/3 [3J, \34 \ . 



But for other values of p, we saw that the 
diffusivity derived from left-passage and de- 
riving function £(£) is incompatible with our 




FIG. 13: (Color online) Main: the variance 
of driving function versus time (which is 
normalized to 1). The slope of line gives 
k ~ 8/3. Inset: the left-passage probability 
versus polar coordinate 0. The solid curve is 
the exact value of P re= g/3 (</>), and the three 
plots correspond to different radii. These 
results are for p = 0. 

previous results in table [H For instance, for 



V 



0.5, we obtain k « 2.24(5 



24]. 



which is 



clearly different from «ew = 4 
Our guess is that the statistics of the span- 
ning curves are affected by the imposed 
boundary conditions in our simulations. To 
see how this can affect our results, we con- 
sidered critical Ising model on a strip geome- 
try (L x /4 = L y ) with two different boundary 
conditions: 

(i) Fixed boundary condition, i.e., all spins 
in the right and left sides of the origin (lo- 
cated at the middle of the lower boundary 
on x axis) are fixed to be up and down, 



11 



respectively. Anti-periodic BC was applied 
along the x direction and the spins at the 
upper boundary were left free. Such BCs 
imposes the interface (domain wall that sep- 
arate spins up from down) emanating from 
the origin to go towards the upper boundary. 
We checked that the spanning interfaces are 
well described by the prediction of SLE with 
k » 3 and d f « 11/8 25, fe 
(ii) Free boundary conditions along both di- 
rections. In this case, the fractal dimension of 
these curves is df ~ 11/8, but the diffusivity 
that is derived from left-passage probability 
and the driving function is k ~ 2. These 
results show that all critical interfaces don't 
satisfy the SLE relation. 

We also simulated critical percolation 



model on a strip geometry with the same two 
boundary conditions. We found the diffusiv- 
ity k ~ 6 37| compatible with fractal dimen- 
sion df ~ 7/4 [38j, for both BCs. This result 
is expected since for the percolation model, 
the statistics of the domain walls are not af- 



fected 
erty 



yy t he BCs due to a locality prop- 



10 



28|. 



In SSM, only for p = 0, we obtain the com- 
patibility between fractal dimension and dif- 
fusivity that is in SAW universality class. It 
is known that the outer perimeter of criti- 
cal percolation clusters in 2c? are described 

n n 

by SAW 39]. In our recent work [121 ]. we 
have shown that the iso-height lines of vari- 
ous discrete models including BD, RSOS and 
EDEN models are in the SAW universality 
class. 
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